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Abstract. We have measured the rapidity distributions dN/dy of 7T , and 
p, p in central Au+Au collisions at ^/s^vjv = 200 GeV. The average rapidity loss 
per participant nucleon is 2.0 ± 0.2 units of rapidity. The strange to non-strange 
meson ratios K/ir are found to track variations of the baryo-chemical potential 
in energy and rapidity. 

In ultra-relativistic heavy ion collisions, final state hadrons are used as a probe 
of the different reaction stages, with a special focus on observables that may reveal 
the existence of an early color deconfined phase, the so-called quark gluon plasma. 
The bulk properties of the collision dynamics are becoming well understood at mid- 
rapidity (\y\ < 0.5). However, in Au+Au collisions at y/spfN = 200 GeV, where the 
beam rapidity is y& = 5.36, this corresponds to less than 10% of the whole rapidity 
range (2 x y^). One of the goals of the BRAHMS experiment [1] is to explore a broader 
rapidity range (—0.1 < y < 3.5). In this paper, we report on some of our latest results 
on identified charged particle yields from the 5% most central Au+Au collisions at 
i/sWn = 200 GeV. 

BRAHMS consists of two 
hadron spectrometers, a mid- 
rapidity arm (MRS) and a 
forward rapidity arm (FS), 
as well as a set of detec- 
tors for global event charac- 
terization [1]. Particle spec- 
tra were obtained by combining 
data from several spectrome- 
ter settings (magnetic field and 
angle), each of which covers 
a portion of the phase-space 
(y,Pr)- The data have been 
corrected for the limited accep- 
tance of the spectrometers us- 
ing a Monte-Carlo calculation 
simulating the geometry and 
tracking of the BRAHMS de- 
tector system. Detector efficiency, multiple scattering and in-flight decay corrections 
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Figure 1. 7r , K and p,p rapidity dis- 
tributions (a) and mean transverse mo- 
mentum (pt) (b). Errors are statistical. 



Charged Hadron Production with BRAHMS 



2 



have been estimated using the same technique, but the data have not been corrected 
for feed-down from resonance and hyperon decays. Details can be found in [1, 2]. 

The pion and kaon spectra were well described at all rapidities by a power law in 
pr, A(l + pT/po)~ n : and an exponential in tot — mo, ^4exp ( mr ~ m ° ), respectively. 
The invariant yields dN/ dy were calculated by integrating the fit functions over the full 
Pt or tot range. The systematic errors on dN/dy, including errors from extrapolation 
and normalization, amount to ~ 10% in the range < y < 1.3 and ~ 15% for y > 1.3. 
Rapidity densities and mean transverse momenta are shown in Fig. 1. Negative and 
positive pions are found in nearly equal amounts within the rapidity range covered. 
In contrast, an excess of K + (p) over K~ (p) is observed to increase with rapidity. 
Figure 1(b) shows the rapidity dependence of (pt)- There is no significant difference 
between positive and negative particles of a given mass. In order to extract full phase 
space densities for tt^ and we have investigated several fit functions: a single 
Gaussian centered at y = (Gl), a sum of two Gaussians (G2) or Woods-Saxon (WS) 
distributions placed symmetrically around y = 0. In Tab. 1 are reported the average 
yields. 



Table 1. Full phase— space yields of -n-^ and extracted from fits to the dN/dy 
distributions (see text). Errors reported here are statistical, systematic errors are 
of the order of 8%. 

7T + 7T~ K + K~ p 

1742 ± 17 1761 ± 16 288 ± 5 241 ± 3 85 ± 4 
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Figure 2. Net— proton ra- 
pidity densities as a function 
of rapidity for three different 
beam energies. 



Figure 3. Net— baryon rapid- 
ity distributions (a) and aver- 
age rapidity loss per partici- 
pant nucleon as a function of 
beam rapidity (b). Errors are 
systematic. 



Figure 2 shows the net-proton densities dN/dy(p) — dN/dy(p) as a function of 
rapidity. Our data [3, 4] (red circles) are plotted together with AGS [5, 6, 7] and 
top SPS [8] data. A drastic change in the shape of the distribution occurs as ^/snn 
increases. At RHIC, the net-proton density is 6.4 ± 0.4 ± 1.0 around mid-rapidity 
(y < 1) and increases to 12.4 ± 0.3 ± 3.2 at y ~ 3. This indicates that collisions 
are much more transparent than at AGS and SPS. The energy available for particle 
production is determined from the energy loss AE of the original nuclei and is related 
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to the degree of stopping. Due to the conservation of the baryon number during the 
reaction, the final state net-baryon (B — B) distribution contains information on AE. 
The nuclear stopping is quantified by the rapidity loss (Sy) , which, for mass symmetric 
collisions, is defined by (Sy) = yb — (y), where (y) is the mean net-baryon rapidity 
after the collision: 



(y) = 



y» dN, 



(B-B) 



(y) 



dy 



(1) 



10 dy 

It is therefore necessary to know the distribution of baryons in the whole rapidity 
range. A Monte-Carlo simulation provided corrections factors for the non-detected 
neutrons, and strange baryons such as A's and S's. Only mid-rapidity hyperon yields 
are known experimentally, so we made assumptions for more forward rapidities. De- 
tails can be found in [3, 4]. Since the data do not cover the entire rapidity range, 
the missing information has been extrapolated from the measured data by fitting the 
latter with a function that is continuous and whose integral is equal to the number of 
participants (357 ± 9 for the 5% most central collisions after a Glauber simulation). 
Figure 3 shows the net-baryon distribution as a function of y (insert) and (Sy) as a 
function yb (main panel). Two functions were used: a sum of Gaussians symmetrically 
positioned around mid-rapidity, and a six order polynomial (dashed lines in insert). 
In both cases, (Sy) amounts to 2.0 ± 0.2. This result is significantly lower than the 
expectation from the empirical scale (Sy) cx yb. The energy loss, calculated after the 
net-baryon fits, is found to be AE — 72 ± 6 GeV per participant, i.e. 25.7 ± 2.1 TeV 
for the corresponding centrality class. 
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Figure 4. Full phase-space K/n ratios as a function of v / sjvjv (a) and rapidity 
systematics at y'sjvlv = 200 GeV (b). The dashed and dotted lines in (b) are 
predictions of the statistical model [9]. Errors are statistical and systematics in 
(a), only statistical in (b). AGS data are from [10, 11], SPS data from [12, 13, 14]. 
Data points at y/sjvjv = 6.3 GeV and 7.6 GeV [12, 13] are preliminary. 

In Fig. 4(a) are shown ratios of the strange to non-strange full phase-space me- 
son yields (K/ir) as a function of ^/sjvat. The ratio K + /ir + shows a fast increase 
from low AGS to low SPS energies (^/snn = 8 GeV), followed by a decrease with in- 
creasing SPS energy (to ^Jsnn — 17 GeV). We find, at ^/snn = 200 GeV, a value of 
0.165 ±0.003 ±0.013, consistent with the SPS ratio at the highest energy. In contrast, 
K~ /n~ increases monotonically but remains below K + /ir + . At ^/s^n — 200 GeV, it 
reaches a value of 0.136 ± 0.002 ± 0.011, which is close to K + /n + at the same energy. 
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The full phase-space ratios are not significantly different from those observed in the 
mid-rapidity region y < 1 (Fig. 4(b)). Indeed, a fit to a straight line in this particular 
range gives K+ /ir+ = 0.156 ± 0.023 (stat + syst) and K~ /w~ = 0.146 ± 0.016. The 
dashed and dotted lines are predictions of the hadron gas statistical model [9] , which 
used a chemical freeze-out temperature T of 177MeV and baryo-chemical potential 
Hb of 29MeV (the authors restricted their fits to yields measured in the rapidity range 
\y\ < 0.5). The agreement with the data is excellent. However, the data deviate from 
the model prediction at higher rapidities, where an increasing excess of K + over K~ is 
observed. This is due to an increase of net-baryon densities [3, 15]. A baryon rich en- 
vironment is favorable for associated strangeness production, e.g. p+p — ► p + K + + A, 
a production channel forbidden to K ~. In the context of the statistical model, this 
translates into an increase of the baryo-chemical potential hb, as already reported 
in [15], where a calculation by Becattini et al. [16] of K + / K~ vs p/p at constant T 
(and varying /is) agrees well with the rapidity dependence of the experimental ratios. 
It is also known that the chemical freeze-out temperature varies strongly in the AGS 
energy range but slightly from SPS to RHIC energies. In this particular energy do- 
main, the K/tt systematics are as well mainly driven by changes in 

In summary, we have measured transverse momentum spectra and inclusive in- 
variant yields of charged hadrons w^, and p,p. The net-proton densities are 
the lowest ever observed in heavy-ion collisions and increase only slightly with in- 
creasing rapidity. The rapidity loss calculated from the net-baryons is found to be 
(Sy) = 2.0 ± 0.2, thus breaking the empirical scaling observed at lower energies. The 
ratio of strange to non strange mesons K/tt ratios are well reproduced by the hadron 
gas statistical model [9] that assumes strangeness equilibration. The increasing differ- 
ence between K + and K~ yields at higher rapidities is explained by a change of the 
baryo-chemical potential \ib with increasing rapidity. 
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